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C/^ . Abstract 

' A short review of the development of the Nambu-Jona-Lasinio (NJL) model is 

CS| ' given. The SU{2) x SU{2) and U{3) x U{3) local quark NJL models are considered. 

. The mechanisms of spontaneous breaking of chiral symmetry and vector dominance 

\^ [ are shown. The local NJL model allows us to describe the mass spectrum and 

' main strong and electroweak decays of the four ground-state nonets of pseudoscalar, 

■ scalar, vector and axial- vector mesons. Applications of this model to the description 

i of mesons in hot and dense medium are discussed. It is shown that for solving 

] problems connected with the description of the meson radial excitations and quark 
confinement it is necessary to consider a nonlocal extension of the NJL model. 

P-i. The main attention is attracted to the description of the methods used in different 
versions of the NJL model. Physical results for low-energy hadron physics obtained 

^ • in these models can be found in the cited works. 



'. 1 Introduction 



The NJL model was proposed in 1961, where the authors attempted to explain 
the origin of the nucleon mass with the help of spontaneous breaking of chiral 
symmetry pp. The model was formulated in terms of nucleons, pions and scalar a 
mesons^. Remind that the fundamental theory of strong interactions, QCD, was 
not constructed at that time. 

After 15 years this model was reformulated at the quark language by the Japanese 
physicists T. Eguchi and K. Kikkawa [HISl. It should be noted that it is supposed 
that all hadrons are formed from constituent quarks with mass m w 300 MeV, 
whereas the QCD theory is based on more light current quarks with mass mo ~ 5 — 7 
MeV. In |S] , it is shown that light current quarks turn into the massive constituent 



^It is worthwhile to note that in the same 1961 two papers devoted to the same problems were published 
in the USSR in Moscow: V.G. Vaks and A.I. Larkin P' and B.A. Arbuzov, A.N. Tavkhelidze and R.N. 
Faustov i3i. 
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quarks due to spontaneous breaking of chiral symmetry. However, they only con- 
sidered the quark NJL model in the chiral limit mo = 0. In this limit all masses of 
pseudoscalar mesons are equal to zero. 

Starting from 1982 M.K. Volkov and D. Ebert with collaborators considered 
a more realistic version of the quark NJL model when uiq 7^ El EI- It al- 
lowed them to describe the mass spectrum, internal properties, and strong and 
electroweak interactions of scalar, pseudoscalar, and vector meson nonets [HI Ell- 
1984, T. Hatsuda and T. Kunihiro applied this model to the description of hadrons 
in hot dense medium jlH I12j . 

After 1986 this model gained popularity. In the following years more than six 
hundred papers devoted to the NJL model were published. Therefore, we cannot 
give here a full list of the corresponding references. We only note some authors who 
gave a noticeable contribution in this field of research. 

The NJL model was especially intensively used for different applications in Ger- 
many (D. Ebert, Humboldt Univ., Berlin; H. Kleinert, Preie Univ., Berlin; H. Rein- 
hardt, Tubingen Univ.; J. Hufner, S. Klevansky, Heidelberg Univ.; W. Weise Mu- 
nich, Tech. Univ.; D. Blaschke, G. Roepke Rostock Univ.; K. Goeke, C.V. Christov 
Ruhr Univ., Bochum; B. Friman Darmstadt GSI and others), Japan (T. Hatsuda, 
K. Yazaki Tokyo Univ.; T. Kunihiro Kyoto Univ.; M. Asakawa, S. Sawada, K. 
Yamawaki Nagoya Univ. and others) and USA (M.D. Scadron Arizona Univ.; 
C.M. Shakin City Coll., N.Y. and others). The NJL model was also studied by 
many physicists from China, England, Prance, Italy, Portugal, South African Re- 
public and other countries. 

In our country this model was actively developed at JINR Dubna (M.K. Volkov, 
G.V. Efimov, M.I. Ivanov, Yu.L. Kalinovsky, A. A. Osipov with collaborators), St.- 
Petersburg State Univ. (A. A. Andrianov and V.A. Andrianov), St. -Petersburg INP 
(D. Diakonov, V.Yu. Petrov) and St .-Petersburg Polytechnic Inst. (A.N. Ivanov). 

It is worthwhile to note that now NJL is widely used for different applications 
as in elementary particle physics as in nuclear physics (see e.g. [T!H ITlj). 



2 SU{2) X SU{2) NJL model 
2.1 Pseudoscalar and scalar mesons 

Firstly, let us consider the simple SU{2) x SU{2) NJL model. The corresponding 
chiral quark Lagrangian is 

G 

^Q, q) = q{x){idx - mo)qix) + {iqix)q{x))'^ + {q{x)iTaj^q{x)f) , (1) 

where q{x) = {u{x),d{x)} are the fields of n, d antiquarks, mo is the current quark 
mass, is the four-quark coupling constant, are the Pauli matrices, and 7^ is 
the Dirac matrix. 

After bosonization the Lagrangian takes the form 

(T^(x) + 7r^(x) 

C{q,q,a,Tr) = q{x){id^ -mo + a{x) + i75TaVra(x))g(x) ——2 , (2) 
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where cr{x), TTa{x) are the scalar and pseudoscalar meson fields. Effective meson 
Lagrangian can be obtained from Lagrangian Q in one quark-loop approximation. 
Here we can see that vacuum expectation value of pion field is equal to zero, whereas 
cr{x) has a nonzero vacuum expectation value (cj)o = ctq. Therefore, it is necessary 
to redefine this field in order a new physical field cr'{x) = cr(x) — ctq have a zero 
vacuum expectation value (cr')o = 0. Excluding linear in cr'(x) terms from the 
Lagrangian we obtain gap equation 



5£ 

6a' 



0, 



mo = m + (To = m(l — SG.,^Ii{m)) 



(3) 



This equation describes spontaneous breaking of chiral symmetry. As a result, light 
current quark mass rriQ turns into massive constituent quark mass m. Note that all 
physical quantities in the NJL model are expressed through the quadratically and 
logarithmically divergent integrals I^[m) and /^(m) 



It{m) 



Nc f 4fc9(A2 _ fe2) ^ ^ 
(27r)4 J m'^ + k^ 47r2 



A2 

A^ - In ( ^ + 1 



(4) 



In 



A2 



+ 1 



1 + 



m 
A2 



These integrals are given in Euclidean space; iVc = 3 is a number of quark colors; 
A is the cut-off parameter. It describes the region where spontaneous breaking of 
chiral symmetry takes place. 

From the Lagrangian the following expressions for the tt and a meson masses 
can be obtained 



2 

dnqq 
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8/f (m) , m: 



ml + 4m^, 



(5) 



where g-j^qq = g^qq = (472^(m))~^/2 ^j-g meson renormalization constants that 
provide a correct coefficient in the kinetic terms for the meson Lagrangian. It is 
easy to see from eqs. (jSJ and (0 that the pion mass is proportional to the first 
power of current quark mass (Gell-Mann-Oakes-Renner relation). As a result the 
pion become the Goldstone particle with zero mass in the chiral limit mg = 0. 
From the weak pion decay tt ^ fiv the Goldberger-Treiman relation follows: 



m 



9nqq 



(6) 



where = 93 MeV is the weak pion decay constant |T 



2.2 Vector and axial- vector mesons 

The quark Lagrangian corresponding to the vector and axial- vector mesons is 



(7) 



3 



After bosonization this Lagrangian takes the form 

C{q, q, p,ai)= g(x)(7^T,p;^(x) + j>^-f'>TaaM^) + ^""^ , (8) 



2G. 



where pa, aia are the fields of vector (p) and axial-vector (oi) mesons. 

Note that in the description of vector and axial-vector mesons a gauge-invariant 
regularization must be used 9^. The quark loop with two vector vertices defines the 
kinetic term of the vector meson and renormalization constant of the vector field 
Qpqq. As a result, the simple relation between g^^gq and gpgg appears [HllHl 

gpqq = V^gaqq (9) 

We obtain for the p-meson mass 



= (10) 



The renormalization constant of the ai meson field Qaiqq coincides with gpgg and 
the mass equals 

M2^=M2 + 6m2 (11) 

2.3 IT — ai transitions 

In the NJL model there are quark loops with pseudoscalar and axial-vector vertices 
that describe n — ai transitions [Qj IT6 | [T71 ITSj . 

These transitions lead to nondiagonal terms in meson Lagrangian of the type 
\/Qmaia{x)dx7ra{x). In order to exclude these terms it is necessary to redefine the 
axial-vector field as 



m . 



ai'ai^) = <a{x) - ^d^T^aix) (12) 

ai 

This leads to additional kinetic terms of the pions and to modification of the renor- 
malization constant g-^qq- Now this constant is not equal to g^qq 

1 - jjrj (13) 

It is interesting to note that allowance for the vr — ai transitions does not affect to 
Goldberger-Treiman relation. 

2.4 Numerical estimations 

Let us define the model parameters. From eqs. ©, ©, fT^ we can define con- 
stituent quark mass thought the observables = 93 MeV, gpgg = 6.14 (this value 
corresponds to the experimental width of p meson) and M^^ = 1.26 GeV pS] 



r- [6 /6 m / / Qm'^\ m 

gpqq = V6g.gg = \l -g^gg = \\ = f \^ - ^^^^ 
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Then 

m^ = ^(^l-^l- ^^p^^ m = 280MeV (15) 

The parameter A can be found from the equation 



9pm = \j^^ ^ A=1.25GeV (16) 

and Gp can be found from equatfons ©, (fTU]) for pion and /o- meson masses. 
Therefore, G^r = 4.9 GeV~'^ and Gp = 16 GeV~^. The value of the current quark 
mass mo is defined from the gap equation ((Hj), niQ = 3 MeV. 



3 /7(3) X /7(3) NJL model 

In order to introduce strange mesons in the model, it is necessary to replace the 
Pauli matrices Tj (i=1..3) by the Gell-Mann matrices Aj (i=0..8, where Aq = ^J^l). 
Let us remind that the Ua{^) problem exists connected with a correct description 
of masses of r/, r]' mesons. Indeed, by using the U{2>) x f/(3) symmetric Lagrangian 
we obtain "ideal" singlet-octet mixing for pseudoscalar isoscalar mesons. Then, one 
of these states contains only u and d quarks and the other state contains only the 
strange quark. This situation contradicts experimental data. 

In order to solve this problem it is necessary to add the t'Hooft interaction JHI 
to NJL Lagrangian j2Ull21l[^l23j . As a result the model for scalar and pseudoscalar 
meson nonets consist of two Lagrangians 

L^JL = ^-(ia-m% + ^^[(g-A,g)2 + (g-i75A,g)2], 

Z:*^ = -i^(det[g(l + 75)'?] + det[g(l-75)g]), (17) 

where q = {■u, d, s} are antiquark fields, mP is a current quark mass matrix with 
diagonal elements m^, m^, (m^ ~ ^S)- 

The Lagrangian (|17j) can be rewritten in the form (see |21j ) 

1 ^ 

C = q{id - m% + - Y.lG'-r^iqWqf + G^+\qi^,>: ^qf] + 

i=l 

+Gi-\qXuq){qXsq) + Gi+\qi-/5Xuq)iqil5>^sq) , (18) 

where 

X'i = X^ (i = l,...,7), A^ = A„ = (V2Ao + A8)/\/3, 
A^ = A, = (-Ao + V2Xs)/V3, 

= G^^^ = G^^^ = G± 4Kmslt{ms), 
Gf ) = g'^^ = Gt^ = gJ^) = G ± AKrUulH^u), 

g1±) = G t 4Km.Ji{ms), G(±) = G, G^t^ = ±4V2Km„/i^(m„). (19) 
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The t'Hooft interaction leads to the additional terms in the gap equations for 
the u, s quark masses 



ruu = m° + SmuGIi (rriu) + 2,2mumsKI^ {mu)Ii (nis) 
rus = ml + SmsGI^{ms) + 2,2K {m.jt{'mu)f 



(20) 



Now let us first consider bosonization of the diagonal parts of the Lagrangian 
H18|) including the isovector and strange mesons. After renormalization of the meson 
fields we obtain 

.2 7 ^2 3 2 7 



£(^,K,ao,i^o*) 



2G^ 



i=l 

—iTr In < 1 



1 



i=4 



2G 



ao 



i=l 



i=4 



id-M 



i=l 



(21) 



/^(+) /^(+) 



I2 {mu,ms) 



Nr 



dtk 



{k^ + m'^){k^ + m^) 



In I — ^ + 1 - m„ In — ^ + 1 
mi / V mf, 



(47r)2(m2 — m^) 

where M is a constituent quark mass matrix, (pi and crj are the pseudoscalar and 
scalar fields. 

As a result, the following expressions for the meson masses are obtained 

1 



M'^ = gl 



Ml = gl 



1 

G^ 

1 



-4[lt{mu) + lHms)] 



+ Z(ms - rriu) 



G„ 



8Ii^(m„) 



+ 4m, 



(22) 



^K* — gK* 



1 



+ (mu + rris 



The nondiagonal part of the Lagrangian has the form 

AC = i {G(+)(g-^75A„(z)' + 2GM(qn5Xuq){qn5Xsq) + G(+)(g-i75A.g)^ 
+Gt)(gA„g)2 +2Gl;)(^-A„5)(g-A,5) + G(-)(gA,g)2} = 

= {qh5^aq)T^i3{qii5^(3q) + {qKq)Tap{q^pq), (" = u, s) (/? = n, , 



+ 



(23) 



where 



"2 Git^ Gf^ 



(24) 



After bosonization we obtain 



-i TrliK^ 1 + 



1 



id-M 

that leads to the following meson Lagrangian 



(25) 



^^"(^ )a/3^/3 7— 



1 



Mi 



9 



(r^)-i-8/f(m„)), 



9vs 



1 



(r^);/ - 8/f (m,) ) , 



(26) 



29vu9vs{T 



P\-i 

us ' 



9a^ [ 



uu 



29an9cTs{T 



S\-l 
us • 



9a^ = 9aqq, 9a s = W2{^s)] 



-1/2 



9riu ~ 9ngqj 9ris 



^1/2 



where au, (Ts are the fields of strange and nonstrange isoscalar mesons, and r]u, ijg 
are the pseudoscalar fields. 

After diagonalization of the Lagrangian (|26)) we find masses of the pseudoscalar 
and scalar mesons tj, tj', a and /o 



M, 



M, 



Mfs + M„^„ T J {MP - MP^f + A{MP,f 



Mt, + M^^ T J(M£ - MS^Y + 4(M„^J2 



(27) 
(28) 



Two additional arbitrary parameters appear in the ?7(3) x U{2>) NJL model: the 
current mass of the strange quark mg and the coupling constant K. We can define 
these parameters using kaon mass and rj - rj' mass difference 



425 MeV, K = 13.3 GeV'^ 



(29) 



Parameters , A and Gt^ remains unchanged, = 280 MeV, A = 1.25 GeV, 
= 4.9 GeV-2 . 
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After that we obtain the following estimations for the masses of the pseudoscalar 
and scalar mesons 



Mr, 



135 MeV, 
520 MeV, 
550 MeV, 



Mk = 495 MeV, 



M 



Ma, = 810 MeV, Mk* 



1000 MeV, 
1130 MeV, 
= 960 MeV, 



(30) 



The experimental data are 

M^o = 134.9764 ± 0.0006 MeV, M^± = 139.6 MeV, 

Mk+ = 493.677 ± 0.016 MeV, M^o = 497.672 ± 0.031 MeV, 



547.45 ± 0.19 MeV, M^, = 957.77 ± 0.14 MeV, 



cro(400-1200) 



400 - 1200 MeV, M 



983.5 ± 0.9 MeV, Mj, 



/o(980) = 
800 MeV. 



980 ± 10 MeV, 



The model parameters are fixed by the masses of pseudoscalar mesons. At the 
same time the scalar meson masses are in a qualitative agreement with experimental 
data. 

Vector and axial- vector mesons in C/(3) x C/(3) version of the NJL model can 
be introduced as in SU{2) x SU{2) version. As a result after bosonization of quark 
Lagrangian we have for vector meson masses |2| 



Mi 



4G. 



M 



K* 



9h 



muf, Ml = (32) 

where Qp = V^QaQ, QK* = V^gK*, 50 = V^gas- Note, the quark loops gives con- 
tribution to mass only for K* meson. As a result the vector meson masses are in 
satisfactory agreement with experimental data 



Mp = 770 MeV, Mk* = 930 MeV, Mj, 



1090 MeV. 



(33) 



It is worthwhile to note that the general Lagrangian for the strong interactions 
of four meson nonets can be expressed in a very compact form 

- 2^ 

12 




M\ 
9 J 



1 2 



My 

7) ^ 



MY Oj, 

g J 2 



where 



Xia\ 



M\ 
9 ) 



n 2 



(34) 



G'^ = d^V" - d^V - ( [V^', V"] _ + [A'',!"] 
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Figure 1: Divergent quark loops with external photons and vector mesons (p. 

The electroweak interactions are introduced in our model in a gauge-invariant 
manner on the basis of the original quark Lagrangian (1). This allows us to describe 
not only the strong processes (strong decays, tttt, ttK scattering and so on) but also 
different electroweak processes such as e.-m. and weak decays, radii, polarizabilities, 
different rare processes(for instance rf — 7r°77). 



4 Vector dominance 

After introducing e.-m. interactions in the NJL Lagrangian the photons can interact 
with the charged mesons only through quark loops. In contrast to the mesons, which 
are composite objects, the kinetic term for the photons is introduced independently 
into the Lagrangian (1). The allowance made for the quark loops leads only to 
renormalization of the electromagnetic fields and the charge. 

The part of the Lagrangian describing the electromagnetic interactions has the 
form 



1- 



id — m 



QA 



(35) 



where 



di/A^, 



(36) 



Q = (^3 -I- Xs/V^)/2 is the operator of the quark charge. 

After calcTilation of the divergent self-energy diagram of the photon (Fig. la), 
we obtain for £em the expression 



id-M 



QA 



(37) 



where 



A' 



1 + 



1^1 
3^ 



1/2 



A' 



1 + 



4e^ 
3^ 



-1/2 



(38) 



Besides the self-energy diagrams involving the photon there are divergent dia- 
grams of mixed type, describing the transitions 7/3^, 70; and 71^. (Fig. lb). The 
inclusion of these diagrams leads to the appearance in the Lagrangian of terms of 
the form 



2* 



1 ^2 , 



fi V I ) 



(39) 
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where pfi„, uj^ y and (\)av tensors constructed from meson fields and derivatives, 
similarly to eq. (|36|) . 

As a result, the part of the Lagrangian that describes the electromagnetic inter- 
actions of the mesons and quarks takes the form 



2 



1 e' ^, 

H F 

2gp 



Ml 
2 
1 



2 

fj. V 



-zTr In 



1 + 



1 



(40) 



where V!* are the fields of vector mesons. 

We diagonalize the kinetic terms by means of the following substitution of the 
fields: 



V2e 



■A' 



(41) 



The electromagnetic field and the charge e' are then renormalized as follows: 

,2\ 1/2 



(42) 



/2' 



1 



4e 
3^ 



-1/2 



/2 



1 



3^ 



1 + 



4e2 
3^ 



-1/2 



It is readily seen that as a result of the two renormalizations [ (jSZI) and (|42|) ] the 
electric charge takes its original value. The final Lagrangian has the form 



1 / ~„2 o -r, ~ 2 

P p!/ + '/'pi/ + ^ pv 



35p, 
zTrln 



1 + 



-^^(7pAv;) 



id-M 2 



(43) 



It is now easy to show that the photons can interact with the charged particles 
only through the neutral vector mesons. In this way, we have automatically obtained 
a model describing vector dominance. Under the sign of the logarithm, the term 
with photons has been completely absorbed by the vector mesons. 

5 Meson in hot and dense matter 



In the last few years the activity in search for a new state of matter, the quark- 
gluon plasma (QGP), has significantly increased. New data are already coming 
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from running experiments on heavy-ion collisions at Brookhaven (RHIC) and CERN 
SPS. New facilities are planned to be constructed to increase our capability in this 
research (LHC, SIS-300). The QGP is expected to reveal itself through modified 
properties of hadronic reactions and their products. 

The NJL model is a very convenient tool for investigation meson behavior in the 
hot and dense matter. First calculations of this type in NJL model were started in 

[in IE]. 

It is possible to use different methods for investigation of meson behavior in the 
hot and dense matter. The most popular one is the Matsubara technique This 
"imaginary time" formalism implies the replacement of the integration over the zero 
component of the momentum by the summation of frequencies 

p° (iuJn + fi) (44) 

where a;„ are the Matsubara frequencies, a;„ = (2n + 1)ttT for fermions and cOn = 
ImrT for bosons; /x is the chemical potential and T is the temperature. 

However, for many applications it is more convenient to use an equivalent rep- 
resentation for the quark propagator derived in the "real time" formalism |29| 



S{p,T,fj.) = ip + m) 



1 

+ 



p2 _ _|_ 



i2Tr6{p'^ - m^){e{p^)n{p,n) + 9{-p^)n{p, -^u))] , (46) 

where 

n{p, ^) = f 1 + exp (47) 



is the Fermi-Dirac function for quarks, E = \/ fP' + m?. This leads to a different 
method for calculation of integrals I^^ {m,T, n), /^^(m, T, /i). First we performed 
contour integration in the complex pQ plane and after that reguralized integrals 
by three-momentum cut-off A3. As a result, divergent integrals I^^ {m^T, fj.) and 
/^^(m,T, /i) take the form 

A3 „ 

A Nc f p 

I-^^{m,T,fi) = j—yjdp—{l-r]{p,^i)-fj{p,^i)), 



A3 

lti^,T,fi) = J dp^{l-iiip,f^)-viP,f^))- (48) 



We defined values of model parameters in vacuum using the same conditions as in 
Sect. 2.4?. Here we assume that model parameters G.^-, Gp, niQ and A3 do not 



2 



Note that different methods of regularization give us some other values of model parameters m = 280 
MeV, A3 = 1.03 GeV, ^ 3.48 GeV"^^ Gp = 15.9 GeV"^ ^nd too = 2 MeV in comparison with sect 
2.4. 
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Figure 2: The behaviour of the quark mass and weak pion decay constant (a), and 
meson masses M^., Mp, Ma^ (b) as function of T. 



depend on T and /i. The dependence of m on T and fi is calculated from the gap 
equation. Then we calculate T and dependence of basic integrals I^^{m, T, /i) and 
/^^(m, T, /i). This allows us to define the dependence on T and /x of all physical 
quantities. 

The behaviour of m{T) is shown in Fig. 2a. For mo = the restoration of chiral 
symmetry is indicated by the vanishing of the order parameter m(T) — > (or quark 
condensate) at critical value of the temperature Tc- When mg ^ the sharp phase 
boundary disappears. Figure 2b exhibits the behaviour of the meson masses Mt^, 
M(j, Mp, Ma^ as function of T. As T increases the mass of the a meson decreases 
sharply with the constituent quark mass. On the other hand, the mass of the pion 
will persistently stay constant until the critical conditions for chiral restoration are 
reached, beyond which it ceases to exist. Mp is merely independent of T, whereas 
shows a sharp decrease similar to that of M^. Above the critical temperature 
one obtains M^- = Mq- and Mp = Ma-^ , as is expected for a chiral symmetric phase. 
In |3U1 131j the value for critical temperature Tc ~ 200 MeV. In lattice recent QCD 
calculations the value for critical temperature Tc ~ 170 MeV . 

Recently, very interesting results have been obtained from the investigation of 
strongly interacting quark matter in the color-superconducting phase. We do not 
consider here the diquark condensation and the related problems. The reader will 
find the details in 

In conclusion, we would like to note that the properties of some particles can 
significantly change when approach to phase transition. The cr-meson, which in 
vacuum is a broad resonance may become very sharp resonance at particular tem- 
perature and chemical potential. When appropriate conditions are reached, this 
can lead to the amplification of some processes mediated by the cj-resonance, such 
as TTvr — > 77, tttt — > vrvr etc. This amplification, if observed in heavy-ion collision 
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experiment, can be interpreted as an event indicating on the approaching the QGP. 



6 First radial excitations of mesons 

In the local version of the NJL model, it is impossible to describe radial excitations of 
mesons. Therefore, for the description of both the ground and first radially excited 
states it is necessary to consider not only the standard local Lagrangian £ but 
also additional nonlocal Lagrangian in Lagrangian £'1°°'°'=^ we introduce a 

form-factor for each quark-antiquark current 

Jj{x)=j j d'^xi(fx2q{xi) Fi{x;xi,X2)q{x2). (49) 

The form-factors -F/(x; xi, 2:2) can be written in a covariant form 35 . Here we 
do not discuss them in details, and only would like to remark the form- factors for 
the ground and first radially excited states can be given in a very simple form in 
momentum space, 

h{k^) = e(A3-|A;^|), (50) 
f2{ki_) = c(l + d|fcx|')e(A3-|fc±|), 

where k is the quark-antiquark pair relative momentum and k± is the part of k 
transversal to the total momentum P 

k± = k- ^P (51) 

The step function, 0(A3 — |A;_|_|), is a covariant generalization of the 3-dimensional 
cut-off in the NJL model. For d < — Aj^^ the form- factor f2{k±) has the form of 
an excited state wave function with a node in the interval < < A3. The 
form- factors H50() are the first terms in a series expansion in k±'^; the inclusion of 
higher radially excited states would require polynomials of a higher degree. The 
factor c describes the change in the strength of the four-quark interaction in the 
radially excited channels relative to the constant Gt^- This constant is defined from 
the mass of radial pion excitation = 1300 MeV. The parameter d is defined from 
the condition^ 

f f d^k /i(k) , , 

// = -iNc / , ■'^^ ' = (53 



•^Here, /„, and I^^ denote the loop integrals with zero, one or two factors f{k±) = f2{k±) in the 
numerator 

- -^^tr / ^/M) (52) 
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Now let us explain the meaning of this condition. Note that in this model we have 
two gap equations 

'J^ = -iN^t, f J^. I (^1)0=0 ,54) 

6ai (27r)4 k - niQ + {ai)o + (a2)o/2(k) 



7a3 (2vr)4 k-mo + ((Ti)o + (^2)0/2(1^) 

In general, the solution of these equations would have (172)0 7^ 0; in this case the 
dynamically generated quark mass, — (cti)o — (cT2)o/i(k) +mo, becomes momentum 
dependent. Condition (|53|) leads to the trivial solution for the second gap equation 
and, as a result, we have only one nontrivial gap equation coinciding with the 
standard gap equation of the local NJL model. 

The free part of the effective action for pions takes the form 

^ = \l^^ jZ<^P)Kt!iP)-)(n (55) 

where i^f/(P) = d^^KijiP) and 

Kii(P) = Zi(p2_M2), K22{P) = Z2{P^ - M^), (56) 



Ki2{P) = K2l{P) = ^/Z^2TP\ 



1 








1 








4/2, 


Z2 = 



with 



Zi = 4/2,Z2 = 4/|^r = ^^/|. 

To determine the physical vr- and vr'-meson states, we have to diagonalize the 
quadratic part of the action. It can be performed with the help of an orthogonal 
tranformation of the fields vri , 112] the details of this procedure can be found in • 
We would like to note that after the expansion in a series of a small current quark 
mass mo (from eq. (|57|) Mf ~ mo), one finds for the physical states 

Ml = Mf + OiMf), (58) 



M. 



2 _ Ml 



i_r2 



i + r^^ + o(M- 



'2 

Thus, in the chiral limit, the effective Lagrangian indeed describes a massless Gold- 
stone boson, the pion tt, and a heavy pseudoscalar meson, vr'. The ratio of the vr and 
vr' weak decay constants can be directly expressed in terms of the physical meson 
masses 

^ - r Ml 
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It is worth noting that the matrix elements for a pseudoscalar meson and the di- 
vergence of the axial-vector current must vanish both for the pion and its radial 
excitation in the chiral limit. In the case of the pion, the matrix element is ~ 
and vanished in the chiral limit as 0. The situation is opposite in the case of 

radially excited state, there F^^' — > in the chiral limit, while M^/ remaining finite. 

Here we discuss only pions. In jSni EH IHH] it is shown that this method can be 
extended to the chiral U{3)xU (3) group for pseudoscalar, scalar and vector mesons. 
In the framework of this model the main strong decays of the radially excited meson 
were described |39| EZj . One of the most interesting results obtained in this model 
concerns the identification of nineteen experimentally observed scalar states for the 
masses 0.4 ~ 1.7 GeV. These states can be interpreted in our model as two scalar 
nonets and a scalar glueball with the mass ~ 1.5 GeV |40| . The first nonet consists 
of the ground state mesons with the masses between 0.4 - 1 GeV. The second nonet 
consists of radially excited scalar mesons with the masses 1.3 - 1.7 GeV. Four scalar 
states and scalar glueball are mixed, as they have the same quantum numbers. 



7 Nonlocal NJL model and quark confinement 

The NJL models have two main defects. They contain ultraviolet (UV) diver- 
gences and do not provide quark confinement. Usually, UV divergences are re- 
moved by using the cut-off parameter A taken at an energy scale of the order of 1 
GeV. The physical meaning of this cut-off is connected with the separation of the 
energy-momentum region, where spontaneous breaking of the chiral symmetry and 
bosonization of quarks takes place. In order to exclude unphysical quark-antiquark 
thresholds from amplitudes of different processes only lowest powers of momentum 
expansion of quark loops are usually used in the NJL models. 

These drawbacks of the standard NJL model can be solved only in the framework 
of nonlocal models. There are many different nonlocal versions of NJL model (see 
e.g. |41[ 1421 1431 1441 145j). Here we demonstrate one version of nonlocal models 
which is motivated by an instanton interaction [IH]- Similar models are considered 

in [HlllHlliniEni- 
The SU{2) X SU{2) symmetric action with the nonlocal four-quark interaction 
has the form 

S{q, q) = j d^x |^(x)(i4 - me)g(x) + ^ (J.(x) J,(x) + J^(x) J^(x)) - 

-^(^^'^(x) Jp^"(x) + > (60) 

The nonlocal quark currents Ji{x) are expressed as 

Jl{x) = j j d!^xid!^X2f{xi)f{x2)q{x-xi)Tiq{x + X2), (61) 

where the nonlocal function f{x) is normalized by /(O) = 1. In (jHOJ the matrices 
Ti are defined as L^ = 1, L^ = i^i^r", T^p" = -f^'r", F^" = 'y^-yf^r". 
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After bosonization the scalar field a will have nonzero vacuum expectation value. 
In order to obtain a physical scalar field with zero vacuum expectation value, it is 
necessary to shift the scalar field. This leads to the appearance of the nonlocal 
quark mass m(p^) instead of the current quark mass mg which can be found from 
the gap equation 

m{p^) = mo + G^--—jf^{p^) d^k - 



= m, + (m(0)-m,)/2(/)^ (62) 

where m(0) is the dimensional parameter which plays the role of the constituent 
quark mass. The quark propagator takes the form 

S{p) = {p - m{p^))-^ . (63) 

We use one of the simplest ansatz for the dynamically generated quark propa- 
gator. In the spirit of |4H I51j we demand that pole singularities are absent in the 
vector part of the quark propagator 

1 l-exp(-pVA') ^g^^ 



"^(P) = ( 1 ■f/x'?\ 7 ) • (65) 



ni?{p'^) + p2 p2 
The expression for m{p) is found to be: 

2 \ 1/2 

exp (p2/A2) — 1 

The mass function m{p^) depends only on one free parameter A, has no any singu- 
larities in the whole real axis and exponentially drops as — > oo in the Euclidean 
domain. From ea. H62() it follows that nonlocal form factors have a similar behavior 
providing the absence of UV divergences in the model. At = the mass function 
is equal to the cut-off parameter A, m(0) = A. From the gap equation we find the 
relation between the four-quark coupling Gi and the nonlocality parameter A 

Moreover, the expression for pion renormalization constant has a simple form 

JVc (3 . C(3) 



where Q is the Riemann zeta function. In the chiral limit one has two arbitrary 
parameters A, G2. We fix their values with the help of the weak pion decay constant 

= 93 MeV, and /o-meson mass Mp = 770 MeV. By using the Goldberger-Treiman 
relation c/^(0) = m(0)/F^ one finds A = m(0) = 340 MeV. 

This simple model leads to reasonable predictions for the a-meson mass M„ = 
420 MeV and strong decay /? ^ vrvr T = 135 MeV. 
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Nonlocal model, in contrast with the local NJL model, can be successfully used 
for the description of not only the constant part of amplitudes of meson interactions 
but also of the momentum dependence of amplitudes at small energies. It is worth 
noting that in the nonlocal model the relative contributions to the pion form-factor 
F^*„+^- and pion radii from the contact diagrams and diagrams with vector mesons 
as intermediate states have reasonable values j52j . In the nonlocal model the contri- 
bution of vector mesons is noticeably smaller than that of the contact diagrams in 
contrast with the local NJL model where these contributions are comparable |53j . 
The vector-meson diagrams play a very important role in the description of the pion 
form- factor Fy^+^- in the time- like region. These diagrams allow us not only to 
describe the /?-meson resonance but also to obtain a correct behavior of the process 
form-factor in the region below 1 GeV. 

8 Conclusion 

Once more emphasize that when the first version of the NJL model was proposed, the 
fundamental theory of strong interactions QCD did not exist. Therefore, that time 
different versions of the phenomenological hadron models were used for description 
of low-energy hadron physics, whereas the description of hadron interactions at large 
energies was very problematic. However, after the construction of QCD and the 
discovery of the phenomenon of asymptotic freedom it became possible to describe 
an interaction of hadrons at large energies by means of perturbation theory. The 
perturbation theory is applicable only for energies larger than 1 GeV when the 
strong coupling constant is smaller than unity. Therefore, for description of the 
low-energy region the usage of different phenomenological models is again needed. 

One of the most attractive models of that kind is the NJL model. The basis 
of the model is the chiral symmetry of strong interactions (as in QCD). The region 
of applicability of this model supplements the QCD perturbation theory. The joint 
application of both the QCD theory and the NJL model allows us to consider the 
whole energy region of the strong hadron interactions. 

M.K. Volkov is grateful to all collaborators (especially to D. Ebert) who took 
part in the construction and development of the NJL model. 
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